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ABSTRACT 

Studies of black hole spin and other parameters as a function of redshift provide 
information about the physical state and merger and accretion histories of the systems. 
One way that black hole spin may be estimated is through observations of extended 
radio sources. These sources, powered by outflows from an AGN, allow the beam power 
and total outflow energy to be studied. In a broad class of models, the beam power 
of the outflow is related to the spin of the black hole. This relationship is used to 
estimate black hole spins for 55 radio sources. The samples studied include 7 FRII 
quasars and 19 FRII radio galaxies with redshifts between 0.056 and 1.79, and 29 
radio sources associated with CD galaxies with redshifts between 0.0035 and 0.291. 
The FRII sources studied have estimated spin values of between about 0.2 and 1; there 
is a range of values at a given redshift, and the values tend to increase with increasing 
redshift. Results obtained for FRII quasars are very similar to those obtained for FRII 
galaxies. A broader range of spin values are obtained for the sample of radio sources 
associated with CD galaxies studied. The fraction of the spin energy extracted per 
outflow event is estimated and ranges from about 0.03 to 0.5 for FRII sources and 
0.002 to about 1 for radio sources associated with CD galaxies; the data are consistent 
with this fraction being independent of redshift though the uncertainties are large. 
The results obtained are consistent with those predicted by numerical simulations 
that track the merger and accretion history of AGN, supporting the idea that, for 
AGN with powerful large-scale outflows, beam power is directly related to black hole 
spin. 
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1 INTRODUCTION 

Black hole spin and mass are the fundamental parame- 
ters that characterize a supermassive black hole. Black hole 
spin is related to the merger and accretion history of the 
black hole as discussed, for example, by Hughes & Bland- 
ford (2003), Volonteri et al. (2005), King & Pringle (2006, 
2007), Volonteri, Sikora, and Lasota (2007), King, Pringle, 
& Hofmann (2008), and Berti & Volonteri (2008). Hughes & 
Blandford (2003) consider the coevolution of black hole mass 
and spin in binary merger scenarios and find that the holes 
are spun down by mergers. Volonteri et al. (2005) study the 
evolution of the full distribution of black hole spins in hier- 
archical galaxy formation theories including prolonged gas 
accretion and binary mergers. Volonteri, Sikora, and Lasota 
(2007) determine and compare spins of black holes in giant 
elliptical galaxies with those of disk galaxies and show that 
elliptical galaxies are likely to have higher spins on average 
than disk galaxies. King & Pringle (2006) suggest that ac- 
cretion onto supermassive black holes occurs in a sequence 
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of randomly oriented accretion episodes. King & Pringle 

(2007) show that black hole growth that occurs in a series of 
small-scale randomly oriented accretion events leads to black 
hole spins that are less than one. King, Pringle, & Hofmann 

(2008) find that accretion of this type leads to supermassive 
black holes with moderate spin, and the spin value of an indi- 
vidual black holes may deviate significantly from the mean 
value. In addition, the mean spin value of the black hole 
population decreases slowly as black hole mass increases. 
Examples of holes with large spin parameter may occur and 
are most likely to be found in giant elliptical galaxies. Berti 
& Volonteri (2008) show that the redshift evolution of black 
hole spins for a given population of sources can be used to 
determine the merger and accretion history of the sources 
and whether the accretion is prolonged or occurs in short- 
lived chaotic events. Thus, studies of spin as a function of 
redshift may be used as a diagnostic of the accretion and 
merger histories of AGN. 

Spins of individual black holes may be determined by 
studies quite close to the AGN, or by studies of outflows 
from the AGN. Each method of measuring black hole spin is 
model-dependent. Spins of three objects have been obtained 
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by studies quite close to the heart of the AGN from X-ray 
measurements. Spins obtained using X-ray measurements of 
AGN are 0.60±0.07 for Fairall 9 (SchmoU et al. 2009), 0.6± 
0.2 for AGN SWIFT J2127.4+5654 (Miniutti et aL 2009), 
and 0.92 - 0.99 for MCG-6-30-15 (Brenneman & Reynolds 
2006; Reynolds & Fabian 2008). The use of extended radio 
sources to study black hole spin is discussed and applied to 
one radio galaxy by McNamara et al. (2009) and to 48 radio 
galaxies by Daly (2009a,b). Daly (2009a) presents a "model- 
independent" method of placing a lower bound on the black 
hole spin. The outflow energy is taken as a lower bound 
on the black hole spin energy, and this bound is combined 
with the black hole mass to obtain a lower bound on the 
black hole spin. Powerful radio galaxies were found to have 
remarkably similar minimum spin values with a weighted 
mean of 0.12 ± 0.01. McNamara et al. (2009) studied the 
radio galaxy MS0735.6-I-7421 in the context of the hybrid 
model of Meier (1999, 2001) and concluded that the outflow 
is likely powered by the spin energy of a maximally spinning 
hole. Daly (2009b) studied 48 extended radio galaxies with 
redshifts between zero and two, including MS0735. 6-1-7421, 
in the context of the Meier (1999, 2001) and Blandford & 
Znajek (1977) models of spin energy extraction and obtained 
an estimate of the spin of each black hole assuming that the 
magnetic field strength is proportional to the black hole spin; 
the black hole spins obtained range from about 0.1 to 1. The 
results suggest that the black hole spin increases slowly with 
redshift over the redshift range from zero to two. 

The work presented here on extended radio sources 
expands and improves upon prior studies, which included 
only radio galaxies, in several ways. A sample of 7 FRII 
quasars is included here, and spin estimates obtained with 
those sources are compared with results obtained with radio 
galaxies. In the prior study, a magnetic field strength B cc j 
was considered. Here, two additional characterizations of the 
magnetic field are considered, and results obtained with dif- 
ferent field strength characterizations are compared. In ad- 
dition, the spin energy per unit black hole mass is obtained 
for the 7 FRII quasars, and a new quantity, the fraction of 
the spin energy extracted per outflow event is obtained and 
studied for 7 FRII quasars, 19 FRII galaxies, and 29 radio 
sources associated with CD galax;ies. 

The methods of estimating the black hole spin, the spin 
energy per unit black hole mass, and the fraction of spin 
energy extracted per outflow event are described in section 
2, results are presented and discussed in section 3, and con- 
clusions follow in section 4. 



2 THE METHOD 

Large-scale radio sources are powered by twin jets that em- 
anate from the AGN. A class of models have been proposed 
and developed in which the jets are powered in part or in 
full by the spin energy associated with a rotating black hole 
and surrounding region (e.g. Blandford & Znajek 1977; Rees 
1984; Begelman, Blandford, & Rees 1984; Punsly & Coro- 
niti 1990; Blandford 1990; Meier 1999, 2001; Koide et al. 
2000; McKinney & Gammie 2004; De Villicrs et al.. 2005; 
Hawley & Krolik 2006). In many of these models there is 
a relationship between the beam power (or energy per unit 
time) Lj carried by the jet, the black hole mass M, spin j, 



and "braking" magnetic field strength B, which takes the 
form Lj (X B^M^j^ where j = a/m, a is defined in terms 
of the spin angular momentum S, the speed of light c, and 
the black hole mass M, a = S/{Mc), and m = GM/(?. 
For example, this proportionality applies to the Blandford 
& Znajek (1977) model (the "BZ" model), the modified BZ 
model discussed by Reynolds, Garofalo, & Begelman (2006), 
and the hybrid model of Meier (1999, 2001). 

Thus, if the beam power and black hole mass for a given 
AGN are known, the black hole spin may be studied for 
various magnetic field strengths (e.g. Blandford 1990, Daly 
2009b): 

3 = n{U4)°-^ Mg-i (1) 

where L44 is the beam power in units of 10'*'' erg/s, B4 is 
the poloidal component of the magnetic field that threads 
the accretion disk and ergosphere in units of lO'' G, and Ms 
is the black hole mass in units of IO^Mq. The constant of 
proportionality k varies by a factor of a few for different 
models; for example, in the hybrid model of Meier (1999) 
K w (1.05)~^^^, while in the model of Blandford & Znajek 
(1977) K»VE. 

A source with no outfiow can still have a significant 
spin, but only sources with outflows provide the information 
necessary to deduce the spin. For example, in the magnetic 
switch model of Meier (1999, 2001), the outflow only occurs 
when the magnetic field strength reaches a particular value. 

The black hole spin j can be used to determine the black 
hole spin energy Es in terms of the black hole mass M (e.g. 
Rees 1984): 

1/2 



Es/Mc = 1 ■ 



1 + 



(2) 



The fraction / of the spin energy extracted during a partic- 
ular outflow event may be empirically determined by taking 
the ratio of the outflow energy to the spin energy: 



/ = E,/Es. 



(3) 



Black hole spins, spin energy per unit black hole mass, 
and the fraction of the spin energy extracted during a par- 
ticular outflow event are estimated here for sources with em- 
pirical determinations of beam power and black hole mass. 
To obtain these estimates, three different characterizations 
of the magnetic field strength are considered. For each char- 
acterization, the range and redshift evolution of black hole 
spin and fraction of spin energy extracted are studied. 

The field strengths considered are an Eddington mag- 
netic field strength, a constant magnetic field strength, and a 
magnetic field strength that is proportional to the black hole 
spin. These three field strengths are related to the black hole 
properties in different ways. Interestingly, the overall results 
and general trends obtained with each characterization are 
rather similar. The Eddington magnetic field strength is the 
field strength such that the energy density of the magnetic 
field is equal to that of a radiation field with an Edding- 
ton luminosity (e.g. Reos 1984; Dermer, Finke, & Menon 
2008); in units of 10"* G, B4,edd ^ 6M~^^'^ . This is the 
field strength that is expected for sources radiating at the 
Eddington luminosity. King (2010) argues that many AGN 
are likely to be radiating at this luminosity for much of their 
lives. To evaluate the impact of an assumed characterizar 
tion of the magnetic field strength on the quantities studied 
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here, we also consider results obtained assuming a constant 
magnetic field strength. This possibility has been suggested 
and considered by several authors (e.g. Rees 1984; Punsly & 
Coroniti 1990; Blandford 1990) and, more recently, Piotro- 
vich et al. (2010). A characteristic value of 10*G, or B4, — 1, 
is generally adopted, and is used here; the results obtained 
here can easily be scaled to any other constant value of the 
field strength. For comparison, a third chaxacterization of 
the magnetic field strength is considered. If it is assumed 
that the fraction of the spin energy extracted per outflow 
event is constant, then observations of FRII radio galaxies 
indicate that the magnetic field strength is proportional to 
the black hole spin, B4 ~ 2.78j (Daly & Gucrra 2002; Daly 
et al. 2009). Interestingly, this is also indicated by a com- 
parison of the beam power predicted in the hybrid model of 
Meier (1999, 2001) and empirical results obtained by Allen 
et al. (2006) and Merloni & Heinz (2007). The results of 
Allen et al. (2006) and Merloni & Heinz (2007) indicate that 
the beam power is proportional to the accretion rate to the 
power 1.3 and 1.6, respectively, while equation (12) of Meier 
(1999) indicates that the beam power is proportional to the 
accretion rate to the power 1.6 when B oc j. Thus, it is pos- 
sible that the field strength is related to the black hole spin, 
and it is useful to consider this case as a comparison. 



3 RESULTS 
3.1 The Samples 

The method is applied to 7 powerful FRII quasars, 19 pow- 
erful FRII galaxies, and 29 central dominant (CD) galaxies, 
and includes most sources for which the black hole mass 
and beam power of the outflow have been empirically de- 
termined. The beam powers and black hole masses for the 
sources are illustrated in Fig. 1 and listed in Tables 1 and 
2. The FRII quasars and galaxies have a similar range of 
beam power and and the quasars have a slightly broader 
range of black hole mass than the galaxies. Most of the CD 
galaxies have substantially lower beam powers and extend 
to lower masses than the FRII sources; this is due to the 
fact that most of the CD sources are quite nearby and can 
be observed to relatively low flux levels. 

The FRII quasars are obtained from the samples of 
Leahy, Muxlow, & Stephens (1989) and Liu, Pooley, & Riley 
(1992), and their properties including beam powers are sum- 
marized by Wellman, Daly, & Wan (1997) and Wan, Daly, & 
Guerra (2000) . All quantities have been computed in a stan- 
dard spatially flat cosmological model with current mean 
mass density fim = 0.3, cosmological constant Q.h = 0.7, 
and a value of Hubble's constant of Ho = 70 km s~^ kpc~^. 
The beam power Lj of the outflow is determined from the 
pressure of the forward region of the radio lobe, the lobe 
width, and the rate of growth of the extended lobes by ap- 
plying the equations of strong shock physics. As noted by 
O'Dea et al. (2009), this determination of the beam power 
does not depend upon the offset of the plasma in the radio 
lobe from minimum energy conditions due to a chance can- 
cellation of the way this offset enters the determination of 
the beam power. The total energy that will be expelled by 
the outflow over its entire lifetime, £*, is obtained from the 
beam power using the relation established for powerful radio 



galaxies E» oc L]''^ (O'Dea et al. 2009; Daly et al. 2009; Daly 
& Guerra 2002). Quantities are obtained for each side of a 
source, and the total value is taken to be twice the weighted 
mean of that for the two sides of the source. 

The sample of 19 FRII radio galaxies is obtained from 
sources studied by Leahy, Muxlow, & Stephens (1989), Liu, 
Pooley, & Riley (1992), Guerra, Daly, & Wan (2000), and 
O'Dea et al. (2009). Their properties including beam powers 
are summarized by O'Dea et al. (2009), and are obtained in 
the same manner as described above for FRII quasars. 

The black hole masses of the FRII quasars are obtained 
from McLure et al. (2006), those for the FRII radio galax- 
ies arc obtained from Tadhunter et al. (2003), McLurc et 
al. (2004), and McLure et al. (2006), and those for the CD 
galaxies are obtained from the values Mbh,Lk listed by Raf- 
ferty et al. (2006) after removing the correction factor of 0.35 
to bring the black hole masses of Cygnus A and M84 into 
agreement with those of Tadhunter et al. (2003) and Mar 
ciejewski & Binney (2001). 

The CD sample is obtained from Rafferty et al. (2006), 
and consists primarily of FRI sources (Fanaroff & Riley 
1974). The total outflow energies Et are estimated by Raf- 
ferty et al. (2006) from the pressure and volume of the cav- 
ity that is occupied by the radio emitting plasma, and the 
beam powers Lj are obtained by dividing the outflow energy 
by the buoyancy timescale t. The authors explain that this 
timescale is probably an overestimate of the true timescale, 
in which case the beam powers would increase, while the en- 
ergy of the cavity could be due to multiple outflow events, 
which would cause the outflow energies and beam powers to 
decrease. 

The source Cygnus A (3C 405) appears in both Table 
1 and 2. The values listed in each table are obtained using 
independent methods. The beam power of (47 ± 8) x lO** 
erg/s listed in Table 1 is obtained by applying the equa- 
tions of strong shock physics to the forward region of the 
source (Wan, Daly, & Guerra 2000; modified to account for 
the different cosmological model adopted), while the value 
of (13 ±7) X 10'"' erg/s listed in Table 2 is obtained by divid- 
ing the energy required for the plasma to create and occupy 
the lobes by the buoyancy timescale (Rafferty et al. 2006); 
Rafferty et al. (2006) explain that the beam power obtained 
in this way may be an underestimate since the buoyancy 
timescale may be an overestimate. Given that there is only 
one source that is common to both methods, it is not pos- 
sible to do a detailed comparison of the methods, and in 
order to compare results obtained with a given method, val- 
ues obtained with that method are not modified. The black 
hole mass estimates for Cygnus A listed in Tables 1 and 2 
are also obtained with two methods, and arc in reasonably 
good agreement. 

The FRII galaxies and quasars studied here are the 
most powerful FRII radio sources at their respective red- 
shift, and the hosts have been identified as massive ellip- 
tical galaxies (e.g. Lilly & Lorigair 1984; Best, Longair, & 
Rottgering 1998; McLure et al. 2004). Thus, the powerful 
FRII and CD sources studied are likely to be drawn from 
the same parent population. Since the FRII sources repre- 
sent the most powerful sources are their respective rcdshifts, 
they likely represent the envelope of the distribution, while 
the CD sample contains sources at low redshift with a very 
broad range of radio powers. 
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Figure 1. Beam powers and black hole masses for 55 AGN; FRII 
quasars are shown as open circles, FRII galaxies are shown as 
solid circles, and CD sources are shown as open stars. 



The radio emitting regions of the sources studied are 
quite large (typically larger than the optical size of the host 
galaxy), and the timescale of the outflows are estimated to 
be millions to hundreds of millions of yeaxs (e.g. Rafferty 
et al. 2006; O'Dea et al. 2009). Thus, it is the very long 
term average properties of the outflows that are studied. 
The beam powers are obtained by applying the strong shock 
equations for FRII sources, and the pressure confinement 
equations for radio sources associated with CD galaxies in 
clusters of galaxies, so beam power determinations do not 
rely on an empirically or theoretically determined relation- 
ship between radio power and beam power. The sources are 
very large, with smoothly varying radio emission; the syn- 
chrotron radiation is not Doppler boosted or beamed, so 
there are no relativistic correction factors of this type. In 
addition, the overall rate of growth of the sources is not rel- 
ativistic (e.g. O'Dea et al 2009; Rafferty et al. 2006). Thus, 
the use of the large extended radio emitting regions around 
AGN to determine the beam power is relatively simple and 
straight forward. In this sense, the AGN outflows studied 
here are quite different from those associated with X-ray bi- 
naries. Outflows from X-ray binaries can be highly variable 
on relatively short timescales, and can be strongly affected 
by Doppler boosting and beaming. The complex nature of 
outflows from X-ray binaries does not allow the direct ap- 
plication of physical principles to obtain the beam power of 
each outflow from detailed studies of that outflow; a theoret- 
ically motivated scaling between radio luminosity and beam 
power is used to obtain the beam power. The different source 
properties studied may make it difficult to compare results 
obtained with large extended radio sources to those obtained 
with X-ray binaries. 



3.2 Determination of the Black Hole Spin 

The black hole spin j depends upon the beam power, black 
hole mass, and magnetic field strength, and is determined for 

each of the three field strengths considered using equation 
(1) in the context of the Meier (1999, 2001) model, labeled 
Jm- Values of Jm and their uncertainties are listed in Tables 
1 and 2, and are shown in Figures 2 - 7. To obtain the frac- 
tional uncertainty of j, the fractional uncertainties of the in- 
put parameters arc combined in quadrature after noting the 
following. When B = Bedd, we have that j oc (L44 /Mg) 
(so the fractional uncertainty of j is rather small). When 
B is constant, the fractional uncertainty of j is larger since 
j oc yTHMg^. When B oc j, j oc {L44y^^Mg;'^^'\ so the 
fractional uncertainty of j is quite small. 

For all three magnetic field strengths, the black hole 
spin increases with increasing redshift. Fits are obtained us- 
ing the FRII galaxies and quasars only. The FRII sources 
are from the complete SCR survey (Bennett 1962) including 
sources from the sample subsets of Laing, Riley, & Longair 
(1983) and Pooley et al. (1987); the selection effects for the 
sample are well understood. As mentioned above, the FRII 
sources studied are the most powerful sources at their re- 
spective redshifts and define the envelope of the distribution 
of sources. 

The slope of Log(jM) as a function of Log(l -I- z) is 
1.1 ± 0.2 for B = Bedd; 0.86 ± 0.36 for B = 10* G; and 
0.43 ± 0.18 for B oc j, and are shown in Figures 2-4. These 
slopes are independent of the value of k and thus are valid for 
a broad range of models. In all cases, the black hole spin is 
increasing with redshift. The reduced x^s of the fits are less 
than one suggesting that the fits may be a bit more signif- 
icant than indicated above. The quasars and radio galaxies 
have a similar range of black hole spin at a given redshift, 
and similar trend with redshift. For the field strengths con- 
sidered, the FRII sources have spins that range from about 
0.2 to 1, and, as noted above, there is a consistent trend of 
spin increasing with increasing redshift. The FRII sources 
are the most powerful sources at their respective redshifts 
and thus are likely to define the envelope of the source distri- 
bution. The black hole spins associated with the CD galaxies 
studied range from about 10~^ to about 0.4. 

Figs. 5, 6, and 7 show the black hole spin as a function 
of black hole mass for the three magnetic field strengths 
considered. The FRII galaxies have a small range of black 
hole mass relative to FRII quasars, and no dependence of 
spin on black hole mass is evident for the radio galaxies. The 
FRII quasars have a broader range of black hole mass than 
the radio galaxies, and the quasar sample is lacking sources 
with masses that overlap those of the galaxies, so it would be 
valuable to acquire and study quasars with black hole masses 
in this range. The CD galaxies have a range of black hole 
mass that is similar to the FRII galaxy and quasar sample, 
and a much broader range of spin, with spins extending to 
very low values. This is likely due to the way that the CD 
sample is selected. 

Note that in the magnetic switch model of Meier (1999, 
2001), outflows are only produced when the magnetic switch 
is activated. Sources with no outflows can have large spin 
values; spin values can only be deduced for sources with 
outflows. 
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Figure 2. Black hole spins obtained with a constant magnetic 
field strength of 10* G; the symbols are as in Fig. 1. The line is 
the best fit to the FRII radio galaxies and quasars only and has 
a slope of 0.86 ± 0.36 and a of 13.28 for 24 degrees of freedom. 
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Figure 3. Black hole spins obtained with an Eddington magnetic 
field strength; the symbols are as in Fig. 1. The line is best fit to 
the FRII radio galaxies and quasars and has a slope of 1.1 ± 0.2 
and a of 18.0 for 24 de grees of freedom. 

3.3 Determination of the Spin Energy and the 
Fraction of Spin Energy Extracted 

After the spin j has been determined, the spin energy Eg 
can be obtained using equation (2), and the fraction of the 
spin energy extracted during a particular outflow event can 
be obtained using equation (3). The spin energy Es depends 
upon the spin j and the black hole mass M, while the frac- 
tion of the spin energy extracted / depends on these quan- 
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Figure 4. Black hole spin obtained with a magnetic field strength 
that is proportional to the spin; the symbols are as in Fig. 1. The 
line is best fit to the FRII radio galaxies and quasars and has a 
slope of 0.43 ± 0.18 and a x^ of 13.2 for 24 degrees of freedom. 




Log (M/10= Mg) 

Figure 5. As in Fig. 2 but for spin as a function of bla<;k hole 
mass. 



titles and the total outflow energy S*. When j is greater 

than unity, Es and / can not be determined (e.g. 3C437, 
3C270.1 and 3C275.1). The values of Es/{Mc^), /, as well 
as Et/{M(?) are are listed in Tables 3 and 4, and / as a func- 
tion of redshift is shown in Figs. 8 and 9. Only two magnetic 
fleld strengths, B = Bedd and Ba = 1, are considered here; 
the conclusion that B oc j was obtained by assuming that 
Et oc Es (or / is constant), as described in section 2, so / 
is, by construction, constant in that case. 

The values of Es/{Mc^) obtained are all much less than 
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Figure 6. As in Fig. 3 but for spin as a function of black iiole 
mass. 
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Figure 7. As in Fig. 4 but for spin as a function of blax;k hole 
mass. 



1. This occurs because r = Es/{Mc^) is a very steep func- 
tion of j, ranging from a value of r w 10~^ for j = 0.1 to 
a maximum value of r w 0.29 for j = 1 (e.g. Rees 1984; 
Blandford 1990), so r is generally small. Equation (2) can 
be rc-arrangod to write the black hole spin j in terms of 
the spin energy per unit black hole mass r = Eg / {Mc?): 



j = 2(2r - 5r^ + 4r^ - r*)^/^. Note that by replacing the 
spin energy Es with the outflow energy E^, this expression 
may used to obtain a lower bound on the black hole spin 
that depends only upon the outflow energy and black hole 
mass (Daly 2009a). 

The fraction of the spin energy extracted per outflow 



Figure 8. Fraction of the spin energy extracted as a function of 
redshift for B = G; the symbols are as in Fig. 1. The best fit 
line to the FRII radio galaxies and quasars is shown and has a 
slope of —0.52 ± 0.41 and a of 14.8 for 21 degrees of freedom. 
This is consistent with no evolution of f(z). 



event / has a smaller range for the FRII sources studied 
than for the CD galaxies studied. The range of / for FRII 

sources is about 0.04 to 0.5 for B = Bedd and 0.03 to 0.3 
for B = 10* G, and FRII quasars have values very similar 
to FRII galaxies. The range of / for CD galaxies is about 
0.02 to greater than 1 for B = Bedd and 0.002 to 0.9 for 
B = lO'' G. The values of / for the CD galaxies may be 
uncertain since this quantity depends primarily upon the 
buoyancy timescale for CD galaxies and, as discussed by 
Rafferty et al. (2006), this is likely to be an upper limit (see 
also the discussion in section 3.1). 

The redshift evolution of / is studied using the FRII 
sources only, and is illustrated in Figs. 8 and 9. For B = 10** 
G, the best fit slope of Log(f) as a function of Log(l + z) 
is —0.52 ± 0.41 with a of 15 for 21 degrees of freedom. 
This is consistent with no evolution of / with z. The values 
of / and redshift evolution of FRII quasars is consistent 
with that of FRII galaxies, though the number of quasars 
is small and more data will be needed to confirm this. For 
B = Bedd, the best fit slope is -1.8 ± 0.4 with a of 276 
for 23 degrees of freedom. The is quite large; to bring the 
reduced x to unity would require increasing the uncertainty 
per point by a factor of about 3.5, which would change the 
uncertainty of the best fit slope to —1.8 ± 1.3, consistent 
with no evolution. Thus, the data are consistent with no 
redshift evolution of / with redshift with a slight hint that / 
may decrease with increasing redshift. As mentioned earlier, 
B oc j was obtained by assuming that / is constant (Daly 
& Guerra 2002; Daly et al. 2009). 
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Figure 9. Fraction of the spin energy extracted as a function of 
redshift for B = Bedd\ the symbols are as in Fig. 1. The best 
fit hne to the FRII radio galaxies and quasars is shown and has 
a slope of —1.78 ± 0.37 and a of 276 for 23 degrees of free- 
dom. To bring the reduced to unity would require increasing 
the uncertainty per point by about 3.5, which would change the 
uncertainty of the best fit slope to —1.8 ± 1.3, so these data do 
not require that /(z) be strongly evolving with redshift. 



4 SUMMARY AND CONCLUSION 

Black hole spins arc estimated for extended radio sources 
fed by outflows from AGN assuming that the beam power 
of the outflow is related to the spin energy of the black 
hole and surrounding region. Samples of 7 powerful radio 
loud quasars, 19 powerful radio galaxies, and 29 CD galax;- 
ies with redshifts between about zero and two are studied. 
All of the sources are likely to be associated with massive 
elliptical galaxies. The black hole spins obtained for power- 
ful FRII quasars and galaxies range from about 0.2 to 1 and 
are largest at high redshift and systematically decrease to 
lower redshift. This is the case for each of the three char- 
acterizations of the magnetic field strength considered. The 
powerful radio galaxies and quasars are drawn from the SCR 
sample and are the most powerful sources at their respec- 
tive redshifts; results obtained with these sources are likely 
to represent the envelope of the distribution. The spins of 
the radio sources associated with CD galaxies range from 
about 0.01 to 0.4. The CD sources are a heterogeneous sam- 
ple and include many nearby low luminosity radio sources. 

The normalization of the black hole spin, k, varies by 
factors of a few for different specific models of spin energy 
extraction from the hole, crgosphere, and surrounding re- 
gion. The trend of black hole spin increasing with redshift 
indicated by the analysis presented here is independent of 
the value of k; difiierent values of k will modify the overall 
normalization of the spins values, changing all of the spins 
by the same factor. If independent spin determinations of 
one or more of these sources are obtained, the value of n 
could be empirically constrained or determined. 

The decrease of the black hole spin with decreasing 



redshift is consistent with the predictions of the models of 
Hughes & Blandford (2003) and King, Pringle, & Hofmann 
(2008). This suggests that beam power is directly related 
to black hole spin. Given the indications of the evolution 
of spin with redshift for massive elliptical galaxies obtained 
here, detailed models for the accretion and merger history 
of the sources and the type of accretion that occurs can 
be constructed, as discussed by Berti & Volonteri (2008). 
Clearly, the results presented here are along the lines of the 
predictions of King, Pringle, & Hofmann (2008) in that the 
mean spin decreases monotonically with decreasing redshift, 
source to source spins have fluctuations of order 0.2, and the 
spin tends to decrease with increasing black hole mass for the 
homogeneous sample of SCR sources (thougfi tfiis is rather 
tentative because the black hole mass is an input to the spin 
determination) . Consistency between model predictions and 
the estimated spins suggests that the beam power is directly 
related to black hole spin. 

The fraction of the spin energy extracted per outflow 
event ranges from about 0.03 to 0.5 for FRII quasars and 
galaxies, and from about 0.002 to about one for radio sources 
associated with CD galaxies. The results are consistent with 
this fraction being independent of redshift. When this frac- 
tion is high, the outflow can signiflcantly modify the spin of 
the hole, while when this fraction is low, the outflow has a 
small effect on the black hole spin. In this case, the merger 
and accretion history of AGN will have a larger impact on 
black hole spin evolution than the outflow history. 

The results obtained here may be considered in light 
of the work of Elvis, Risaliti, & Zamorani (2002), who find 
that black holes associated with AGN that produce the X- 
ray background must be rapidly spinning, on average, at 
the time that background is produced. Models developed 
by Gilli, Salvati, & Hasinger (2001) suggest that the X-ray 
background is produced by sources with a range of redshift 
with much of the background produced by sources with red- 
shift between 1 and 3. The elliptical galaxies studied here are 
found to have high spins at redshift of about 2, with spins 
slowly decreasing to lower redshift. If the sources that pro- 
duce the X-ray background have spins that evolve in a man- 
ner similar to those studied here, the results of Elvis, Risaliti, 
& Zamorani (2002) might suggest that the most mild evo- 
lution of spin with redshift is favored; in the context of the 
current study, this would favor the magnetic field strength 
B oc j, which yields spins that evolve as (l-|-2)'^ *^^°'^*. Note 
that this field strength is also favored by a comparison of the 
empirical results of Allen et al. (2006) and Merloni & Heinz 
(2007) with spin energy extraction models, as described at 
the end of section 2. Another factor that would affect the 
production of the X-ray background is the normalization of 
the spin, through the factor k. If the sources that produce 
the X-ray background have spins similar to those studied 
here, it would suggest a value of k ~ \/5 as indicated by the 
BZ model. In this case, the spins of FRII sources vary from 
values of order unity at a redshift of about 2 to about 0.7 
at a redshift of about (e.g. Daly 2009b). Given the range 
of redshift of sources contributing to the X-ray background, 
and the indications from this study that there is a range of 
spin values at a given redshift for radio sources, a detailed 
study would have to be carried out to determine whether 
the sources that produce the X-ray background could have 
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spin normalizations and redshift evolution similar to those 
obtained here for radio sources. 

The spin values obtained here for AGN are in good 
agreement with typical spin values reported for X-ray bi- 
nary systems, which arc summarized by Fender ct al. (2010). 
Studies of X-ray binaries have not indicated a link between 
black hole spin and the properties of the jets (Steiner et 
al. 2010; Fender et al. 2010; Gou et al. 2010), except for 
a correlation between the transient jet normalization and 
spin (Fender et al. 2010). According to Fender et al. (2010), 
the lack of a relationship between jet properties and spin 
for X-ray binaries could be due to a problem with the jet 
properties or speed determinations; a problem with the spin 
measurements; or a real lack of dependence of jet properties 
on spin. Studies of AGN suggest a link between the beam 
power of large-scale outflows and black hole spin (Sikora, 
Stawarz, & Lasota 2007). To date, the samples studied are 
rather small, and very different parts of the systems arc 
used to study the jet properties of X-ray binaries and the 
beam powers of AGN, as described in section 3.1. Thus, 
it is not clear how significant the differences are between 
the relationship of beam power and spin for AGN and the 
relationship of jet properties and spin for X-ray binaries. 
If differences truly exist, it might suggest that more than 
one mechanism may lead to the formation of jets from the 
vicinity of a black hole. In related work, Fernandes et al. 
(2010) studied powerful radio galaxies, concluding that the 
accreted energy is channeled into jets with maximum effi- 
ciency for the most powerful sources, and suggesting that 
black hole spin plays a dominant role in the production of 
jets. Punsly (2010) compared jet efficiencies with currently 
popular "no net flux" models of spin energy extraction and 
found that these models could not account for the empiri- 
cally determined jet efficiencies. The detailed study of Ko- 
vacs, Gergely, & Biermann (2010) suggests that efficient jets 
can be produced in the context of the BZ model of spin en- 
ergy extraction. 

Overall, the results obtained with powerful 3CR quasars 
are consistent with those obtained with SCR galaxies. The 
sample of quasars studied is small and larger samples of 
quasars arc needed to quantify this result. The quasar sam- 
ple studied has a black hole mass range that is larger than 
that of the SCR galaxies studied and a dearth of black hole 
masses that overlap those of the galaxies. There is a hint 
that the quasars may lie at the edges of the distribution of 
spins for SCR sources at a given redshift, and this might be 
related to black hole mass. This can be tested with larger 
quasars samples, and with samples of powerful radio galaxies 
with a broader range of black hole mass. 
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Table 1. Black Hole Spins of FRll Sources 
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RG 
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31 ±8 


14± 10 


1.6 ±0.5 


0.25 ± 0.09 


0.38 ± 0.26 


0.37 ±0.13 
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RG 


0.63 


20 ±6 


9.1 ±6.2 


2.0 ±0.7 


0.25 ±0.09 


0.48 ± 0.34 


0.41 ±0.14 


3C34 


RG 


0.69 


65 ±9 


16± 11 


1.5 ±0.5 


0.33 ±0.12 


0.48 ± 0.34 


0.41 ±0.15 


3C441 


RG 


0.707 


65± 12 


18± 12 


1.4 ±0.5 


0.32 ±0.12 


0.45 ± 0.32 


0.40 ±0.14 


3C 55 


RG 


0.72 


180 ± 50 


14± 10 


1.6 ±0.6 


0.58 ±0.22 


0.91 ±0.66 


0.57 ±0.21 


3C 247 


RG 


0.749 


35 ± 9 


26± 18 


1.2 ±0.4 


0.19 ±0.07 


0.22 ±0.16 


0.28 ±0.10 


3C 289 


RG 


0.967 


85 ± 19 


27 ±21 


1.2 ±0.4 


0.30 ±0.12 


0.34 ±0.26 


0.35 ±0.13 


3C 280 


RG 


0.996 


53± 15 


27 ±21 


1.2 ±0.4 


0.23 ±0.1 


0.26 ± 0.21 


0.31 ±0.12 


3C 356 


RG 


1.079 


250 ± 85 


28 ±22 


1.1 ±0.5 


0.50 ± 0.22 
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^ Total beam powers arc the weighted sum of the beam power from each of the two lobes of each source. 
Input values for each side of each radio galaxy (RG) are obtained from by O'Dea et al. (2009), who used the 
data sets of Leahy, Muxlow, & Stephens (1989), Liu, Pooley, & Riley (1992), Guerra, Daly, & Wan (2000), 
cUid new observations. The beam power is obtained by applying the equations of strong shock physics to the 
source. Input values for radio loud quasars (RLQ) are obtained from the compilations of Wellman, Daly, & 
Wan (1997) ajid Wan, Daly, & Guerra (2000) after converting quantities to the cosmological model adopted 
here; the data used for the study originated in the observations of Leahy, Muxlow, & Stephens (1989) and 
Liu, Pooley, & Riley (1992). Only sources with black hole mass estimates are included. 

The value for 3C 405 (Cygnus A) is obtained from Tadhunter et al. (2003); values for the next four sources 
are obtained from McLure et al. (2004); values for the remaining 14 RG are obtained from McLure et ai. 
(2006). Values for the RLQ are obtained from McLure et al. (2006). 
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Table 2. Black Hole Spins of Soiircos in Galaxy Chistors 
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J- 1 ^ y 


1 4 lb 4 


All _|_ A AC 
U. 1 -L ^ U.UO 


1 « -L n 10 
u. 10 m u. -LU 


94 -1- 07 

U.Z'X ^ U.U 1 


A133 


CD 


0.06 


6.2 ± 1.4 


20± 10 


1.3 lb 0.3 


0.093 ± 0.025 


0.12 ±0.06 


0.21 ±0.05 




ou 


U.UOo 


1 -U 1 /I 
l.D it 1.4 


OQ _l_ 1 1 
Zo it 11 


1.0 it U.o 


A A/l /I _l_ A AOO 
U.U44 It U.UZZ 


U.U04 ± U.UoO 


A 1 /I -U A At; 
U. 14 it U.Uo 


A2029 


CD 


0.077 


0.87 ±0.27 


60 ± 36 


0.77 ±0.23 


0.020 ±0.007 


0.015 ± 0.009 


0.073 ± 0.023 






n 0S1 


1.0 ±0.5 


26 ± 14 


1.2 ±0.3 


A A'^'i -1-0 019 
U.UOo It U.UIZ 


A AQQ -1- A nO'i 
U.UOo III U.UZO 


A 1 9 _|_ A Ozl 
U. IZ It U.U4 


A2597 


CD 


0.085 


0.67 ±0.58 


8.6 ±2.9 


2.0 ±0.3 


0.047 ± 0.022 


0.093 ± 0.051 


0.18 ±0.05 


3C 388 


CD 


0.092 


2.0 ± 1.8 


17 ±7 


1.4 ±0.3 


0.057 ±0.028 


0.081 ± 0.049 


0.17 ±0.05 


PKS 0745-191 


CD 


0.103 


17 ±9 


31 ±16 


1.1 ±0.3 


0.12 ±0.04 


0.13 ±0.07 


0.21 ±0.06 


Hercules A 


CD 


0.154 


3.1 ±2.5 


20± 11 


1.3 ±0.4 


0.066 ± 0.032 


0.086 ± 0.060 


0.17 ±0.06 


Zw 2701 


CD 


0.214 


60 ±62 


17 ±9 


1.4 ±0.4 


0.31 ±0.18 


0.44 ± 0.32 


0.40 ±0.14 


MS 0735.6+7421 


CD 


0.216 


69 ±51 


20± 11 


1.3 ±0.4 


0.31 ±0.14 


0.41 ± 0.28 


0.38 ±0.13 


4C 55.16 


CD 


0.242 


4.2 ± 3.0 


14 ±7 


1.6 ±0.4 


0.09 ±0.039 


0.14 ±0.09 


0.22 ±0.07 


A1835 


CD 


0.253 


18± 13 


54 ±36 


0.81 ±0.27 


0.10 ±0.05 


0.076 ± 0.057 


0.16 ±0.06 


Zw 3146 


CD 


0.291 


58 ±42 


74 ±53 


0.70 ± 0.25 


0.15 ±0.07 


0.10 ±0.08 


0.19 ±0.08 


Almost all of the 


sources 


have FRI 


or amorphous radio structure except for 


a few exceptions 


such as 





Cygnus A (Birzan et al. 2008). 

Beam powers arc obtained from RaiTerty et al. (2006) who combine the total cavity pressure and volume 
with the buoyancy timescale to determine the beam power. 

Black hole masses arc obtained from the values Mbh.Lj^ listed in Table 3 of Raffcrty ct al. (2006) after 
removing the 0.35 correction factor introduced in that paper. This brings the black hole mass estimates 
of Cygnus A (3C 405) and M84 into reasonably good agreement with the independent determinations of 
Tadhunter et al. (2003) and Maciejewski & Binney (2001), respectively. 



R. A. Daly 



Table 3. FRll Black Hole Spin Energy and Fraction of Spin Energy Extracted 



Source 


type 




Es{b=Bedd) 


f(^B — ^edd) 


Ba(S = 10' 


'G) 


/(B = lO-^G) 


Mc'-' 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


3C 405 


RG 


1.3 ±0.4 


0.0067 ± 0.0022 


0.195 ±0.016 


0.0092 ± 0. 


,0055 


0.141 ±0.041 


3C 244.1 


RG 


1.9 ± 1.4 


0.0050 ± 0.0037 


0.379 ± 0.050 


0.019 ± 0. 


.027 


0.102 ±0.072 


3C 172 


RG 


3.6 ±2.6 


0.014 ±0.011 


0.253 ± 0.033 


0.07 ±0. 


.13 


0.049 ± 0.034 


3C 330 


RG 


3.3 ±2.4 


0.023 ± 0.018 


0.145 ±0.019 


0.07 ±0. 


,12 


0.049 ± 0.034 


3C 427.1 


RG 


1.8 ± 1.3 


0.0077 ± 0.0057 


0.235 ± 0.029 


0.019 ±0, 


,027 


0.096 ±0.067 


3C 337 


RG 


2.5 ± 1.8 


0.0079 ± 0.0059 


0.316 ± 0.043 


0.032 ± 0, 


,048 


0.079 ±0.055 


3C34 


RG 


2.4 ± 1.7 


0.014 ±0.011 


0.168 ±0.012 


0.031 ± 0, 


,048 


0.077 ±0.054 


3C441 


RG 


2.3 ± 1.7 


0.013 ±0.010 


0.173 ±0.015 


0.027 ±0, 


,041 


0.086 ±0.061 


3C 55 


RG 


4.3 ±3.2 


0.048 ± 0.042 


0.089 ±0.012 


0.16 ±0, 


,44 


0.027 ±0.020 


3C 247 


RG 


1.2 ±0.9 


0.0047 ± 0.0037 


0.253 ± 0.034 


0.0063 ± 0, 


,0094 


0.19 ±0.14 


3C 289 


RG 


1.6 ± 1.3 


0.0113 ± 0.0094 


0.141 ±0.016 


0.015 ±0. 


,024 


0.109 ± 0.085 


3C 280 


RG 


1.3 ± 1.1 


0.0070 ± 0.0059 


0.186 ± 0.026 


0.009 ± 0. 


,014 


0.15 ±0.12 


3C 356 


RG 


2.7 ±2.3 


0.034 ± 0.033 


0.079 ± 0.013 


0.042 ± 0. 


,078 


0.064 ± 0.052 


3C 267 


RG 


2.6 ±2.2 


0.029 ± 0.027 


0.089 ± 0.012 


0.042 ± 0. 


,078 


0.062 ± 0.051 


3C 324 


RG 


1.6 ± 1.4 


0.015 ±0.014 


0.105 ±0.019 


0.014 ±0, 


,024 


0.116 ±0.099 


3C 437 


RG 


5.5 ± 5.1 


0.16 ±0.27 


0.0350 ± 0.0045 








3C 68.2 


RG 


2 ± 1.9 


0.022 ± 0.023 


0.091 ± 0.015 


0.022 ± 0, 


,043 


0.091 ±0.085 


3C 322 


RG 


3.3 ±3.2 


0.065 ± 0.077 


0.0510 ± 0.0072 


0.07 ±0, 


,16 


0.047 ± 0.045 


3C 239 


RG 


2.8 ±2.8 


0.052 ± 0.062 


0.054 ± 0.010 


0.05 ±0, 


,11 


0.060 ± 0.059 


3C 270.1 


RLQ 


13.1 ±7.7 












3C 268.4 


RLQ 


2.1 ± 1.3 


0.052 ± 0.064 


0.041 ± 0.010 


0.027 ± 0. 


,055 


0.080 ± 0.050 


3C 68.1 


RLQ 


1.33 ± 0.78 


0.019 ±0.019 


0.0710 ± 0.0090 


0.008 ± 0. 


,015 


0.167 ± 0.098 


3C 175 


RLQ 


0.71 ± 0.41 


0.0057 ± 0.0055 


0.125 ±0.012 


0.0025 ± 0. 


,0046 


0.29 ±0.17 


3C 254 


RLQ 


1.9 ± 1.1 


0.011 ±0.012 


0.164 ±0.025 


0.020 ± 0, 


,040 


0.094 ± 0.056 


3C 275.1 


RLQ 


7.7 ±4.6 


0.016 ±0.016 


0.488 ±0.071 








3C 334 


RLQ 


0.73 ± 0.43 


0.0044 ± 0.0042 


0.168 ±0.020 


0.0030 ± 0, 


,0056 


0.25 ±0.15 



The total outflow energy Et is taken to be the weighted sum of the outflow energy of each side of a source; 
the values for the RG were obtained from Daly (2009a) who used the outflow energies listed by O'Dea et al. 
(2009) for the radio galaxies; values for RLQ were obtained in an identical manner. 
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Table 4. CD Bl;«:k Hole Spin Eiioigy and Fraction of Sjiin EnergA' Extracted 



Source 




Mc'-' 


f{B — Bedd) 


£;s(S=10*G) 


f(B = lO^G) 




(2) 


(3) 




(4) 




(5) 


(6) 


M84 


0.00019 ± 0.00023 


0.000010 + 0, 


.000011 


0.019 + 0, 


,022 


0.00010 + 0.00012 


0.0020 + 0.0024 


M87 


0.00052 ± 0.00028 


0.000024 + 0, 


,000013 


0.0210 + 0, 


,0089 


0.000097 ± 0.000077 


0.0050 + 0.0017 


Centaurus 


0.0016 ± 0.0010 


0.000030 + 0, 


,000018 


0.052 + 0, 


,029 


0.00012 + 0.00010 


0.0130 + 0.0035 


HCG 62 


0.0018 ± 0.0022 


0.000024 + 0, 


,000028 


0.075 + 0, 


,082 


0.00014 + 0.00021 


0.0130 + 0.0059 


A262 


0.0034 ± 0.0020 


0.000039 ± 0, 


,000024 


0.086 + 0, 


,012 


0.00016 ± 0.00013 


0.0210 + 0.0076 


Perseus 


0.25 + 0.19 


0.00030 ± 0, 


,00018 


0.81 + 0, 


,53 


0.00061 ± 0.00057 


0.41 + 0.32 


PKS 1404-267 


0.0047 + 0.0045 


0.00012 + 0, 


,00012 


0.038 + 0, 


,010 


0.00073 + 0.00097 


0.0060 + 0.0034 


A2199 


0.083 + 0.057 


0.00047 + 0, 


,00034 


0.178 + 0, 


,035 


0.00081 + 0.00083 


0.103 + 0.049 


A2052 


0.022 + 0.022 


0.00030 + 0, 


,00030 


0.072 + 0, 


,013 


0.00061 + 0.00075 


0.036 + 0.016 


2A 0335+096 


0.017 + 0.013 


0.000058 + 0, 


,000046 


0.293 + 0, 


,037 


0.00014 + 0.00016 


0.122 + 0.063 


MKW 3S 


0.99 + 0.65 


0.0017 + 0. 


,0011 


0.591 ± 0, 


,005 


0.0068 ± 0.0060 


0.146 + 0.049 


A4059 


0.023 + 0.018 


0.00012 + 0, 


,00010 


0.200 + 0, 


,062 


0.00014 + 0.00017 


0.167 + 0.098 


Hydra A 


1.24 + 0.74 


0.00131 + 0, 


,00062 


0.95 + 0, 


,44 


0.0040 + 0.0032 


0.314 + 0.035 


A85 


0.0093 + 0.0078 


0.000045 + 0, 


,000037 


0.21 + 0, 


,14 


0.000054 + 0.000064 


0.173 + 0.082 


Cygnus A 


0.65 + 0.46 


0016 + 


,0011 


0.412 + 0, 


,010 


0.0019 + 0.0021 


0.34 + 0.17 


Sersic 159/03 


0.32 + 0.27 


0.0016 + 0, 


,0014 


0.204 + 0, 


,027 


0.0032 + 0.0039 


0.102 + 0.053 


A133 


0.27 + 0.15 


0.00108 + 0, 


,00059 


0.247 + 0, 


,062 


0.0019 + 0.0019 


0.144 + 0.070 


A1795 


0.046 + 0.046 


0.00024 + 0, 


,00025 


0.188 + 0, 


,013 


0.00037 + 0.00049 


0.125 + 0.063 


A2029 


0.018 + 0.012 


0.000050 + 0, 


,000034 


0.353 + 0, 


,031 


0.000029 + 0.000035 


0.62 + 0.37 


A478 


0.013 + 0.010 


0.00014 + 0, 


,00010 


0.096 + 0, 


,048 


0.00018 ± 0.00022 


0.072 + 0.040 


A2597 


0.093 + 0.085 


0.00027 + 0, 


,00025 


0.344 + 0, 


,061 


0.0011 + 0.0012 


0.086 + 0.033 


3C 388 


0.067 + 0.068 


0.00041 + 0, 


,00040 


0.17 + 0, 


,15 


0.0008 + 0.0010 


0.083 + 0.030 


PKS 0745-191 


0.49 + 0.34 


0.0019 + 0, 


,0013 


0.259 + 0, 


,038 


0.0021 ± 0.0023 


0.24 + 0.12 


Hercules A 


0.34 + 0.34 


0.00054 + 0, 


,00053 


0.639 + 0, 


,010 


0.0009 + 0.0013 


0.37 + 0.21 


Zw 2701 


4.5 + 5.2 


0.013 + 0, 


,015 


0.36 + 0, 


,38 


0.026 + 0.040 


0.175 + 0.084 


MS 0735.6+7421 


18 + 16 


0.012 + 0, 


,012 


1.44 + 0, 


,25 


0.022 + 0.031 


0.82 + 0.49 


4C 55.16 


0.19 + 0.16 


0.00102 + 0, 


,00089 


0.182 + 0, 


,047 


0.0025 + 0.0030 


0.076 + 0.043 


A1835 


0.19 + 0.19 


0.0012 + 0, 


,0011 


0.17 + 0, 


,12 


0.0007 + 0.0011 


0.26 + 0.17 


Zw 3146 


1.1 + 1.2 


0.0027 + 0, 


,0028 


0.42 + 0, 


,31 


0.0013 + 0.0020 


0.91 + 0.61 



The total energy is obtained from Railerty et al. (2006). 



